Broadband polarization retarder 
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We derive a new kind of effective broadband polarization device using the analogy between the 
equation for the polarization of light, propagating through an optically anisotropic medium, and the 
optical Bloch equation in adiabatic regime, which describes coherent laser excitation of a two-state 
atom. 
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Two simple and very convenient tools to describe the 
change of the polarization of light, transmitted through 
anisotropic optical media, are the Jones vector and 
Stokes vector [l|-|3( • The equation of motion for the Jones 
vector in a medium with zero polarization-dependent loss 
has a Schrodinger equation form d, [j| , while the equa- 
tion of motion for the Stokes vector has a Bloch equation 
form These properties have been used recently to 

draw analogies between the motion of the Jones vector, 
alternatively the Stokes vector, and a spin-1/2 particle in 
nuclear magnetic resonance, or alternatively an optically 
driven two-state atom in quantum optics, res pec tively 
described by Schrodinger or Bloch equations |9|-|l4| . 



The phase delay between the fast and the slow eigenpo- 
larizations of traditional birefringent polarization devices 
scales in proportion to the frequency of the light and thus 
such devices are frequency dependent: a half wave plate 
and a quarter wave plate work for a single frequency. 
At the same time, there exists a broadband polarization 
transforming device, known as Fresnel rhomb [U-IH], but 
the quality of polarization transformation is limited by 
the dispersive properties of the material, the rhomb is 
made of. Except the Fresnel rhomb, broadband polar- 
ization converters have proven to be inefficient, whereas 
efficient polarization converters are narrowband. 



The adiabatic polarization conversion proposed here 
is frequency independent: any input polarization state 
is transformed to the desired output polarization state 
regardless of the wavelength. This conversion acts as a 
broadband device, which is limited only by the absorptive 
characteristics of the device instead of its birefringence 
bandwidth. It is noteworthy that the proposed technique 
is analogous to the adiabatic technique in quantum optics 
fl5l - [l9l |. and hence enjoys the same advantages in terms 
of efficiency and robustness. 



We begin with the derivation of the basic equations of 
light polarization evolution in anisotropic nonmagnetic 
linear crystal with negligible absorption. For such media, 
where no free charges and no currents are present, the 



Maxwell equations are 
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where E is the electric vector, B is the magnetic vector, 
e is the electric permittivity tensor and c is the light 
velocity in vacuum. Writing the Maxwell equations (p} 
in a Cartesian coordinate system for which z direction, 
the direction of the propagation of light, coincides with 
one of the optical axes of the crystal, then the tensor e 
can be written as [TH, [l4| 
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where we assume that it is a symmetric function of z only. 
The components e xx (z) , e yy (z) and e xy (z) are expressed 
in terms of the principal values e x (z) and e y (z): 
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where (p is an angle of rotation between the xy plane of 
the chosen Cartesian coordinate system and the optical 
axes. 

If we take the curl of Eq. (fTc)) , using the vector identity 
V x (V x E) = V (V • E) -AE, and if we assume that 
the electric field varies slowly in the plane transverse to 
the propagation direction, V • E w 0, we arrive at the 
wave equation for the electric field: 



V 2 E 



<9 2 E 



(4) 



We shall deal with light propagation in one dimension. 
Let a monochromatic plane wave with frequency uj propa- 
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gate through the inhomogeneous crystal along the z axis: 

" A x (z) e ( ikz - w€ . " 
A y (Y) e^ lk 



E(z,t) = 
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Since the electric field E (z, t) depends only on the longi- 
tudinal coordinate z, we can replace V 2 by d 2 /dz 2 , and 
as a result we obtain the following simplified form of Eq. 
(|4"|), written as two scalar equations: 
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It is generally acceptable to neglect the second term on 
the left-hand side of each equation, as this term is much 
smaller compared to the third one. This is known as the 
slowly varying amplitude approximation poj and is valid 
whenever we have 

d 2 A x 



dz 2 
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This condition requires the fractional change at a dis- 
tance of the order of an optical wavelength to be much 
smaller than unity. In this case Eq. (|6]) becomes 
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where we have taken into account that k 2 = (oj 2 e z ) /c 2 . 
We note that the repeating diagonal term e z can be re- 
moved from the last equation by incorporating it as an 
identical phase in the amplitudes A x and A y . 

Considering Eq. ([3]), we rewrite the last equation as 
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If we map coordinate dependance into time dependance, 
Eq. @ is equivalent to the Schrodinger equation for 
a two-state atom in rotating-wave approximation 
Il9j . where A x and A y are the probability amplitudes for 
the ground state (horizontal polarization) and the excited 
state (vertical polarization). The off-diagonal element CI 
in Eq. ([9]) is known as Rabi frequency, while the element 
A corresponds to the atom-laser detuning [TJI, [l6| . 

We can express the (real-valued) components of the 
Stokes vector 1 3] by the components of Jones vector: 
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Since the medium has negligible absorption, the value 
of Sq is conserved and therefore Eq. ^ turns into the 
optical Bloch equation [TJ|, [l6| : 
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Traditional polarization devices work on resonance 
(A = 0), which is achieved by adjusting the plane of the 
incident light, so that it makes an angle ip = n/ 4 with 
the fast axis (see Eq. (jlObjl ). Furthermore, to achieve 
good polarization conversion an exact area condition is 
required (j Cldz = mr with n being half integer). Con- 
sequently, the traditional retarders depend on the prop- 
agation length (the thickness of the plate) and are not 
broadband. Hence they are not optimal. To this end, we 
propose an alternative and robust adiabatic technique to 
create broadband polarization retarder. 

Let us assume that we use layers of the same material 
to make our retarder and that we change only the angle 
ip from layer to layer. Then we can write Eq. (1121) . using 
Eq. (fTU|) . in the so-called adiabatic basis 
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The connection between the Stokes vector S (z) in the 
original basis and the Stokes vector S A (z) in the adia- 
batic basis is given by 

S(z) = R(z)S A (z), (14) 

with R(z) being the unitary transformation matrix 



sin(2iy9) v / 2cos(2(/3) sin(2yj) 
cos (2ip) — v2sin(2(/?) cos(2ip) 
i -i 



(15) 



For adiabatic evolution of the system there are no tran- 
sitions between the amplitudes S A , S A and S A . Hence 
\S A \ , \S A \ and remain constant p^4l9j . Mathe- 

matically, adiabatic evolution means that in Eq. f| 13[) 
the non-diagonal terms can be neglected compared to 
the diagonal terms, which holds when we have [15l4l9j 
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Thus, adiabatic evolution requires smooth z-dependance 
of the angle tp and large rotary power For a pure adi- 
abatic evolution, the solution of Eq. (|13p is very simple, 



S A {z } ) = U A {z fl z l )S A {z l ) 

where the adiabatic propagator U A {zj,z 
and contains only phase factors: 



V A {z h z l ) = 
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with the adiabatic phase 77 = f* f \idz. The propagator or explicitly: 
in the original basis, the reduced Mueller matrix, is 

U(z f , Zi ) = R(z f )U A (z f ,z i )R J <{z i ), (19) 



cos (a) cos (/3) + sin (a) sin (/3) cos (77) — sin (a) cos (/3) + cos (a) sin (/3) cos (77) sin (/3) sin (77) 
U(zf,Zi)= cos (/3) sin (a) cos (77) — cos (a) sin (/3) cos (a) cos (,0) cos (77) + sin (a) sin (/3) cos {(3) sin (77) , (20) 



sin (a) sin (77) 



cos (a) sin (77) 



cos (77) 



with a = 2ip (zi) and (3 = 2tp (2/). 



Equation (|20|) gives the general adiabatic evolution sce- 
nario for an arbitrary polarization. Naturally and with- 
out loss of generality we can fix x and y axis of the Carte- 
sian coordinate system to coincide with the fast and slow 
optical axis for the first layer of the crystal, which is 
equivalently to set tp(zi) = 0. In this case Eq. (|2"0"|) re- 
duces to simpler form: 



U(z f ,; 



cos (/3) sin (/3) cos (77) sin ((3) sin (77) 
— sin ((3) cos (/3) cos (77) cos (/3) sin (77) 
— sin (77) cos (77) 



(21) 

Obviously there are two terms in the last equation that 
do not depend from the adiabatic phase 77, therefore those 
terms are not frequency dependent and a retarder, that 
use only those two terms, will be frequency independent. 

Broadband rotation for linearly polarized light: 
If initially the light is linearly polarized in horizontal di- 
rection, S(zi) — (1,0,0), then from Eq. (|2"Tj) we have 



Si{z f ) = cos {2ip (z f )), 
S 2 {z f ) = - sin (2<p (z f )) . 
S 3 (z f ) = 0. 



(22) 
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Therefore we can end up with a vertically polarized light, 
S(zf) = (—1,0, 0), if we set the final angle <p(zf) = n/2. 
This process is reversible: if we start with a vertically 



polarized light and set the final angle to be tp{ z f) = 7r /2, 
we achieve reversal of the direction of motion and we end 
up with a horizontal polarization. 

Analogously, if initially the light is linearly polarized in 
horizontal direction, S(z^) = (1, 0, 0), but now we choose 
the final angle to be <p(zf) = 3ir/4 or f{zj) = ir/4, 
then we end up either with a linear +45° polarized light, 
S(zf) — (0,1,0), or with a linear —45° polarized light, 
S(zf) = (0, —1,0). Again, those processes are reversible. 

The presented adiabatic retarder is advantageous be- 
cause the polarization rotation depends only on the an- 
gle of rotation. These retarder is frequency independent 
and it is robust against variations of the propagation 
length, rotary power, etc., in contrast to the traditional 
retarders. 

In conclusion, we have shown that using the anal- 
ogy between the equation, which describes the polar- 
ization state of light, propagating through an optically 
anisotropic medium, and the optical Bloch equation in 
adiabatic regime, which describes coherent laser excita- 
tion of a two-state atom, it is possible to build effective 
and broadband polarization retarder. 
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